Violation of the entropic area law for Fermions 
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We investigate the scaling of the entanglement entropy in an infinite translational invariant 
Fermionic system of any spatial dimension. The states under consideration are ground states and 
excitations of tight-binding Hamiltonians with arbitrary interactions. We show that the entropy of 
a finite region typically scales with the area of the surface times a logarithmic correction. Thus, 
in contrast to analogous Bosonic systems, the entropic area law is violated for Fermions. The re- 
lation between the entanglement entropy and the structure of the Fermi surface is discussed, and 
it is proven, that the presented scaling law holds whenever the Fermi surface is finite. This is in 
particular true for all ground states of Hamiltonians with finite range interactions. 



Entanglement is a phenomenon of common interest in 
the fields of quantum information and condensed mat- 
ter theory. It is an essential resource for quantum in- 
formation processing and intimately connected with ex- 
citing quantum phenomena like superconductivity, the 
fractional quantum Hall effect or quantum phase transi- 
tions. Crucial to all these effects are quantum correla- 
tions, i.e., the entanglement properties, of ground states. 
These have recently attracted a lot of attention, leading 
to new insight into quantum phase transitions and renor- 
malization group transformations and triggeringthe 
development of new powerful numerical algorithms [2j ■ 

A fundamental question in this field is concerned with 
the scaling of the entropy — which is for pure states syn- 
onymous with the entanglement. That is, given a ground 
state of a translational invariant system, how does the 
entropy of a subsystem grow with the size of the consid- 
ered region? Originally, this question appeared first in 
the context of black holes, where it is known that the 
Bekenstein entropy Q is proportional to the area of the 
horizon, which led to the famous conjecture now known 
as the holographic principle 0, Q . The renewed inter- 
est, however, comes more from the investigation of spin 
systems and quantum phase transitions. Moreover, the 
scaling of the entropy is of particular interest concerning 
the choice of the right ansatz-states in simulation algo- 
rithms. 

In the last years, especially one-dimensional spin 
chains have been studied extensively and it is now be- 
lieved that the entropy diverges logarithmically with the 
size of a block if the system is critical, and that it sat- 
urates at a finite value otherwise For a number of 
models 0, @, in particular those related to confor- 
mal field theories in 1+1 dimensions [ïfl lïïj. this could 
be shown analytically revealing a remarkable connection 
between the entropy growth and the universality class of 
the underlying theory. At the same time the diverging 
number of relevant degrees of freedom provides a sim- 
ple understanding of the failure of DMRG methods for 
critical spin-chains. 

For several spatial dimensions a suggested entropic 
area law ^3] could recently be proven for the case 



of a lattice of quantum harmònic oscillators (quasi-free 
Bosons), where again the entropy grows asymptotically 
proportional to the surface. On heurístic grounds this 
can be understood from the fact the system is non- 
critical: an energy gap gives rise to a finite correlation 
length, which in turn defines the scale on which modes 
inside the subsystem are correlated with the exterior. Al- 
though a general area law for gapped lattice systems has 
not been proven so far, the case of quasi-free Bosons is 
often considered as paradigmàtic. In fact, recently de- 
veloped simulation algorithms based on ansatz-states ex- 
hibiting the presumed entropy scaling are highly promis- 
ing 0. 

The fact that in some 1-d systems a vanishing energy 
gap leading to a diverging correlation length results in 
the logarithmically diverging entanglement entropy in- 
evitably raises the question about the behavior of gapless 
systems in more than one dimension. 

The present paper is devoted to the study of the entan- 
glement entropy in gapless Fermionic systems of arbitrary 
spatial dimensions. We establish a relation between the 
structure of the Fermi sea and the scaling of the entropy 
and prové that a finite non-zero Fermi surface implics 
that the entanglement grows proportional to the surface 
of the subsystem times a logarithmic correction, i.e., 



S~L d -HogL 



(D 



if the system under consideration is a c?-dimensional cube 
with edge length L. Thus, in contrast to analogous 
Bosonic systems the entropic area law is violated for 
Fermions. 

Before we start to prové this result a brief discussion 
of the notion of locality — necessary for the concept of 
entanglement — is in order. In spin systems as well as in 
the Bosonic case of harmònic oscillators the tensor prod- 
uct structure of the underlying Hilbert space naturally 
leads to an unambiguous notion of locality. In the ab- 
sence of such a tensor product structure either in general 
quantum field theory settings 01 or m the present case 
of Fermions |lfi| , one has to identify commuting subalge- 
bras of observables and assign them to different parties. 
In our case the relevant àlgebra is the one spanned by 
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the Fermionic creation and annihilation operators cj and 
Cj satisfying the usual anti-commutation relations. We 
assign the modes with j = 1, . . . , n to one party A (the 
interior) and the other modes j = n + 1, n + 2, . . . to the 
other party B (the exterior). Then parity conservation or 
even stronger, particle conservation 0, leads to a super- 
sclcction rule, which implies that all physical operators 
acting on A commute with those acting on B, leading to 
well-defmed notions of locality and entanglemcnt. 

Let us now introduce the prerequisits for the proof. 
Consider a number preserving quadratic Hamiltonian 



H 



E 



,P C k C 



T = T \ 



(2) 
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describing Fermkms on a d-dimensional cúbic lattice, so 
that each component of the vector indices a, corre- 
sponds to one spatial dimension. Translation symmetry 
is reflected by the fact that T is a Toeplitz operator, i.e., 
T a ,/3 = T a -p depends only on the distance between two 
lattice points. The Hamiltonian J2J is diagonalized by a 
Fourier transform leading to the dispersion relation 



e(k) =J2 T « 



k e [~tt,tt} c 



(3) 



All thermal and excited states of the Hamiltonian H are 
Fermionic Gaussian state [ÏÜj , which are completely char- 
acterized by their correlation matrix 



7c*/3 



5 a p - 2tr [pc^C/j] 



(4) 



The correlation matrix describes a pure state iff j 2 = 1, 
so that all eigenvalues of 7 are ±1, and the ground state 
correlation matrix is given by 7 = tÇt. Ground states for 
different Fermion densities are then obtained by adding 
a chemical potential, i.e., replacing T by T + fj,l. If we 
characterize the Fermi sea by the corresponding indicator 
function 6(k) £ { 0, li , then the respective correlation 
matrix is given by [20j 



1 



7a/3 



(2tt) 



dki 



dk d [l-29{k)} e lk - {a ~ í3) 



(5) 

Note that this characterizes not only ground states but 
pure Gaussian states in their most general form (as long 
as they obey particle conservation and translation sym- 
metry) . 

The state of a subsystem, e.g., a cube with edge length 
L, is described by the corresponding L d x L d sub-matrix 
of 7, which we will denote by 7. This subsystem can be 
decomposed into normal modes by a canonical transfor- 
mation from U{L d ) such that the state of each normal 
mode has a Fock space representation of the form 



1 - A, 



|1>(1| 



1 + A, 



(6) 



where the A^ are the eigenvalues of 7. The entropy of the 
subsystem can then be expressed as 

L d 

S(j) = £ h(Xj) , (7) 

Kx) = — log — — log . (8) 

Since a direct computation of S(j) via the diagonaliza- 
tion of 7 is yet highly non-trivial in the simplest one- 
dimensional case with nearest neighbor interaction 
one relies in general on finding good bounds on the en- 
tropy. We will use quadratic bounds on h(x) of the form 
f(x) = o(l - x 2 ) + b 21]. The best lower bound is given 
by a = 1, 6 = leading to 



S(7) >tr [I-7 2 ] 



(9) 



The set of tight quadratic upper bounds can be parame- 
terized by the point Xq £ [0, 1) for which /(xo) = h(xo) 
become tangent We will couple this bound to the 

block size L via xq = 1 — l/g(L), where g(L) = L/logL. 
Straight forward but lengthy calculations show then that 
the entropy as a function of L is asymptotically upper 
bounded j23j by 



5(7)<0(tr [I-7 2 ] logff(£)) • 



(10) 



Hence, together with the lower bound this means that 
tr [l — 7 2 ] essentially determines the asymptotic sealing 
of the entropy. 

The necessity of coupling the upper bound to L can 
easily be understood physically, when one recalis that 
there is always a choice of the local basis in which each 
normal mode inside the block is only correlated with at 
most one mode outside |24j. With increasing block size, 
more and more modes inside lose their correlations with 
the exterior, such that the number of nearly pure normal 
modes dominates more and more. Since the correspond- 
ing eigenvalues are Xj ~ ±1, the point Xq, for which the 
bound is tight, should tend to 1 as L — > 00. Setting 
xq = 1 right from the beginning is, however, not possible 
since the derivative of h(x) at this point diverges. 

Let us now investigate the sealing of tr [l — ■y 2 ] . To 
this end we introduce the positive Fejér kernel 



Fl(x) 



cos(Lx) — 1 
cos(a;) — 1 



(11) 



and we will abbreviate F L {h) by F L {k). Then fol- 
lowing Eq.© we have 

= j^ïdj dkdk' 8(k)[l-6(k')]F L (k-k') 



= (2^2d Jdq^(q)F L (q) , 
E(q) = j dk 9(k)[l - 9(q + k)] . 



(12) 
(13) 
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To further evaluate Ea. 1112(1 we have to exploit the fact 
that with increasing L the Fejér kernel Fl(x) becomes 
morè and more concentrated around x = 0. In fact, 
F L (0) = L 2 , f*^dxF L (x) = 2ttL and for all e > there 
cxists a finitc constant c f such that 



í dqZ(q)F L (q)< c e + í dq E(q) F L 

J[-TÏ,TT] d J[-£,e] d 



(q) • (14) 



The crucial point here is that c e does not depend on 
L. Hence, the asymptotic scaling of the entropy de- 
pends only on the behavior of the function S(ç) in an 
e-neighborhood of the origin. 

The function S(q) has a very intuitive interpretation: 
it is the volume of the part of the Fermi sea in fc-space, 
which is no longer covered if we shift the Fermi surface 
by a vector q (see Fig.l). So what is the behavior of 
S(g) near the origin? Obviously, H(0) = and E(q) > 0. 
Moreover, S will typically not be differentiable at q = 0, 
but it will rather have the structure of a pointed cone. 

Let us assume that the Fermi sea is a set of non-zero 
measure with a finite non-zero surface. This means in 
particular, that almost all points with 9(k) = 1 are in- 
terior points of the Fermi sea and it implies that S is a 
continuous function in fc-space j2^|. In fact, an infinite 
boundary could lead to a discontinuity of S at the origin. 
This restriction excludes both, trivial cases (zero entropy 
due to zero surface) and exòtic cases (fractal or Cantor 
set like Fermi seas). For all other cases it enables us to 
bound E(q) in a neighborhood of the origin by pointed 
cones in the following manner: consider the surface of 
the closed interior of the Fermi sea in one unit celi of 
the reciprocal lattice. Let s(q) be the area of the pro- 
jection of this surface onto the hyper-plane with normal 
vector q, where we account for each front of the sur- 
face. If, for example, the Fermi sea consist out of two 
disjoint three-dimensional sphcres with radius r, then 
s(q) = 2nr 2 in every direction. Since S(g) is the vol- 
ume in which the Fermi sea changes upon shifting it by 
q, we have in an e-neighborhood of the origin that E(q) is 
given by s(q)| |q| I2- Using the fact that the Fermi surface 
is assumed to be finite we know that s(q) is boundcd 
from above by a finite constant s + . Let us assume for 
the moment that there exists a non-zero lower bound 
s~ as well. Then we can bound the integral in Eq. i|14|) 
by replacing E(q) with s ± ||q||2- Exploiting further that 
Iklli > IMI2 > ||ç||i/Vo and that Fl is symmetric, leads 
to upper and lower bounds which are up to a finite con- 
stant given by 

2 d s ± í dqF^WqW, = 2 d s ± V f dqF L {q) qi (15) 

J[0,7v] d ~[ J[0,iv] d 

= 2ds ± (2irL) d - 1 í dxF L (x)x . 
Jo 

The remaining integral is the Fejér sum of a linear func- 




FIG. 1: Consider the Fermi surface and shift it by a vector 
q in fc-space. H(ç) is then given by the (dark gray) area of 
the Fermi sea which is no longer covered by its translation. 
The scaling of the entanglement entropy depends only on the 
behavior of E in the vicinity of q — 0. 



tion, which can be evaluated [23j to 

dxF L (x)x = 2(l + c 7 +hL2 + 7/>(£) + 0(ir 1 )) 

= 21ni + 0(l) , (16) 

where c 7 ~ 0.577 is Euler's constant and ip denotes the 
digamma function. 

We still have to discuss the case inf ç s(q) = 0. In this 
case we have to use different linear bounds for different 
directions in Ea. (|15fl . Since s(q) will be larger than some 
s~ > at least in one direction, we are in the end led to 
the same type of integral and thus to the same asymptotic 
scaling. 

Putting it all together we have indeed that tr [l — 7 2 ] 
scales as L^ 1 ^ lnL since all the involved constants are 
finite and depend not on L but merely on the structure 
of the Fermi sea. The above argumentation holds undcr 
the assumption that the Fermi surface is not too exòtic. 
However, if the interactions are finite in range, then the 
Fermi surface of the ground state is differentiable to infi- 
nite order and it is in particular finite. In general, how- 
ever, one has to check whether or not the structure of the 
Fermi sea gives rise to an infinite slope or a discontinuity 
of E(q) at the origin. 

The existence of Fermi surfaces leading to a scaling of 
the entanglement entropy which surpasses the above law, 
can easily be understood: consider a Fermi sea given by 
a checkerboard with squares of edge length l. Then shift- 
ing the Fermi surface by l along any lattice axis yields 
5 = 2ir 2 such that a naive limit l — * would indeed 
give rise to a L d scaling of the entropy. Needless to 
say that the checkerboard does not have a well-defined 
limit — however, following the same idea, more sophisti- 
cated Cantor set like constructions will do the same job 
without any caveat. In fact, for d = 1 such states were 
constructed in |2l·|. 

Remarkably, fractal or Cantor set like structures are 
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known to appear in tight binding models. The most 
prominent example is the Azbel-Hofstadter Hamiltonian 
[26| with non-integer flux, leading to the famous Hofs- 
tadter butterfly for the spectrum. Since the interaction 
matrix Tki — expi f. A(s)ds (with A being the vector po- 
tential) is quasi-periodic and not translational invariant, 
this case is, however, not directly covered by the above 
argumentation. The question, which physically interest- 
ing translational invariant Hamiltonians give rise to a 
violation of the above scaling law via a fractional Fermi 
sea, remains an interesting problem for future research. 

In conclusion we derived a method of relating the struc- 
ture of the Fermi sea in tight-binding models to the scal- 
ing of the entanglement entropy. For every finite non-zero 
Fermi surface we proved the violation of a strict area law 
(as it is assumed for non-critical systems 0,0,113) by 
a logarithmic correction, i.e., 

c^L^HogL < S < c+L 4 - 1 (log Lf , (17) 

with constants c± depending only on the Fermi sea. By 
the strong sub-additivity of the entropy the same scaling 
behavior holds true also for other regions, e.g. spheres, as 
long as they can be nested into two cubes of edge lengths 
L and cL with c independent of L. The additional log_L 
in the upper bound is presumably an artefact (cf. [jUl^) 
coming from the incompatibility of tight quadratic upper 
bounds with the binary entropy function at ±1. 

Note finally, that the derived result can be applied 
to spin models in one dimension 0, . In this case a 
Jordan-Wigner transformation maps Fermionic operators 
onto Pauli spin operators such that every tight-binding 
Hamiltonian with nearest neighbor interactions in Eq. J2J) 
is then mapped onto a spin Hamiltonian of the form 

OC 

with some couplings hj,. Conversely, every such Hamilto- 
nian is covered by Eq.©, and we are in general allowed 
to add arbitrary interaction terms differing from those 
in Eas. (|18ll9() by a sequence of <r z s in between every 
two Pauli operators. For higher dimensions, however, an 
analogous construction fails, since then Jordan-Wigner 
transformations do no longer preserve locality. 
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